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Abstract 

The energy and material processing industries are traditionally characterized by very large- 
scale physical capital that is custom-built with long lead times and long lifetimes. However, 
recent technological advancement in low-cost automation has made possible the parallel op¬ 
eration of large numbers of small-scale and modular production units. Amenable to mass- 
production, these units can be more rapidly deployed but they are also likely to have a much 
quicker turnover. Such a paradigm shift motivates the analysis of the combined effect of lead 
time and lifetime on infrastructure investment decisions. In order to value the underlying real 
option, we introduce an optimal multiple stopping approach that accounts for operational 
flexibility, delay induced by lead time, and multiple (finite/inhnite) future investment oppor¬ 
tunities. We provide an analytical characterization of the hrm’s valne function and optimal 
stopping rule. This leads us to develop an iterative nnmerical scheme, and examine how the 
investment decisions depend on lead time and lifetime, as well as other parameters. Further¬ 
more, our model can be used to analyze the critical investment cost that makes small-scale 
(short lead time, short lifetime) alternatives competitive with traditional large-scale infras¬ 
tructure. 
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1. Introduction 

The energy and material (e.g. water and petrochemicals) processing industries are tra¬ 
ditionally characterized by very large unit-scale physical capital. For example, individual 
electric power generators rated in the 100s of MW, distillation towers in refineries measuring 
100,000s of barrels-per-day of capacity, and mining trucks capable of hauling 400 tons of ore 
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are common sizes in their respective industries. However, the recent emergence of low-cost 
automation technologies makes a modular, small-scale, and potentially distributed approach 
possible with comparable aggregate production capacity. This calls for a re-examination of the 
current “bigger-is-better” paradigm. Indeed, abandoning economies of unit scale in favor of 


econo mies of mass-production presents several opportunities, as discussed in iDahlgren et al. 


(120131) . In contrast to the very long lifetimes of 25 years or more for typical large-scale capital, 
small and mass-produced equipment might be endowed with a much shorter physical lifespan. 
This could be either by design in construction or in operation with a limited maintenance 
schedule. Furthermore, with a smaller unit scale comes the ability to build to stock and 
drastically shorten lead times between the investment and operation. These factors provide 
the hrm with additional flexibility to engage and disengage a given activity, which should be 
accounted for in the investment valuation and decision. 

Motivated by this anticipated paradigm shift, we introduce here a framework that incor¬ 
porates operational flexibility, lead time, capital lifespan, and multiple (hnite/inhnite) future 
investment opportunities. To this end, we formulate an optimal multiple stopping problem, 
where the hrm maximizes the expected discounted reward from sequential investments. In 
particular, the project’s reward function captures the operational hexibility of temporarily 
suspending production to avoid negative cash hows. Additionally, it depends explicitly on 
two crucial elements, namely, lifetime and lead time. Under capacity constraint, the hrm’s 
consecutive investments are separated (or refracted) by the capital lifetime. Hence, the hrm’s 
investment decision bears similarity to the valuation of a forward-starting swing call option 
written on the reward. 

Our main result is the characterization of the hrm’s value function and optimal stopping 
rule. The optimal timing for multiple investments is described by a sequence of critical price 
thresholds. Furthermore, these thresholds are shown to be decreasing and converge to that 
corresponding to the case with inhnite investment opportunities. Our analysis lends itself to 
an iterative algorithm that numerically solves for the optimal value function and all exercise 
thresholds, as oppos ed to the simulation approach commonly found in existing li t eratu re for 
swing-type options ( Meinshausen and Hamblv . 2004 Chiara et ah . 2007 : Bender . 2011 ). We 
examine the impacts of lead time and lifetime, as well as other parameters, on the investment 
decisions. Moreover, our model is also useful for analyzing the critical investment cost that 
makes small-scale (short lead time, short lifetime) alternatives competitive with traditional 
large-scale infrastructure. 

As is well known, the real option approach to investment decisions can increase the 
project value above and beyond those rooted in classi cal net present value (NPV) argu¬ 
ments ( Sick and Gamba . 2010 : Dixit and Pindv^ . 1994 ). In its inherent myopia, the NPV 
approach falls short of capturing the flexibility of timing since it, at best, only gives the deci¬ 
sion maker an indication of whether or not a single investment should be made at the current 
time. Standard real investment option analysis addresses this limitation. Examples of the 
implementatio n of such an analysis r e gardin g individual investme n t deci si ons can be found in, 
among others. iKaslow and PindvckI (119941): iFraver and Uluderel (120011) : iLumlev and Zervos 
( 20011 ): Carelli et al. ( 2010l) anZ westner and Madlenei ( 2012 ). Nevertheless, the valuation of 
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multiple consecutive investments has not received a great deal of attention in the industries 
in question. In the current paradigm, the individual investment is typically both large and 
long-lived, so the hrm’s next investment decision, including capacity replacement, will be an 
issue of the far future. With such a long horizon, discounting would reduce future cash flows 
to minimal present values. However, with shorter lead times and lifetimes, future investments 
can potentially have signihcant bearing on current capital budgeting decisions, as we will show 
in this paper. 

Our valuation approach accounts for multiple and possibly inhnite future investment op¬ 
tions, which in turn allows for a comparison of capital investments with different lifetimes. 
To illustrate, let us compare two projects with different lifetimes, namely, one with 10 years 
and another with 3 years. To address the difference in investment horizon, one approach 
would be to consider a longer horizon of 30 years, the smallest common multiple of 3 and 10, 
and compare the corresponding net present values. However, this method implicitly assumes 
seamless consecutive investments at the end of each lifetime. This may not be optimal to the 
hrm since it ignores the embedded timing option in each investment decision. 

Another feature of our valuation framework is to incorporate the flexibility to delay future 
investments depending on market conditions. Another advantage over the pairwise compar¬ 
ison is that we can easily examine how sensitive the investment timing is with respect to 
different model parameters, especially for all lead times and lifetimes. We do acknowledge 
that durability of capital is a path-dependent variable, leading to a variable lifetime. In con¬ 
trast, our approach uses a hxed lifetime independent of the frequency of operation, and it 
may favor the posited status quo of capital with a longer prescribed lifetime. Our choice of 
a hxed lifetime, albeit a simplihcation, results from the trade-off between model tractability 
and realistic details. 

The use of swing options has been most prevalent in energy delivery contracts where the 
holde r has the right to alter, or ‘swing’ , volum es up or down at the start of each time pe¬ 
riod f Jaillet et ah . 2004 : Deng and Oren . 2006h . Rather than hxed-period contracts, the only 
constraint on the timing in our setting is a minimum refraction time between consecutive 
investments. As such, the valuation of the swing contract can be forrnulated as a refracted op¬ 
timal multiple stopping problem (see e.g. ICarmona and Touzil (120081) : ICarmona and Davanik 
f 2008l H. Other hnancial applicat i ons involving multiple exercises inclu de employee stock op¬ 
tion valuation (iLeung and Sircarl. l2009l: iGrasselli and Hendersonl . 120091 ). and the operation of 
a physical asset flLudkovskil . 120081) . 

This paper is structured as follows. In Section [2] we introduce the investment decision 
as a general optimal multiple stopping problem. We also introduce a reward function that 
incorporates the flexibility to temporarily shut down production to avoid a negative cash flow. 
In the following section we state and prove sufficient conditions on a general reward function 
for the multiple stopping problem to have a well dehned stopping boundary. We also present 
an algorithmic approach to hnding the solution. In Section 0] we further motivate and analyze 
the specihc reward function in the context of basic infrastructure investments and present 
numerical results. These results include a comparison of investment scenarios with short 
lifetimes and short lead times versus a traditional scenario with comparatively long lifetimes 
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and lead times. Finally, conclnding remarks and potential extensions are discnssed in Section 

El 


2. Problem Formulation 


We consider a firm that has the ability to invest in capital eqnipment that prodnces a single 
good in a given market. Fnrthermore, we assnme that this hrm is acting as a price taker in 
this market, with a hnite capacity constraint. These conditions imply that the investment 
decision will be of ‘bang-bang’ type, that is, either invest np to this limit or not at all. From 
an operational cost perspective, we are agnostic to whether the total capacity comes in one big 
nnit or in 10,000 smaller nnits. Under snch circnmstances, we can withont loss of generality 
analyze the investment timing on a per-nnit-capacity basis. 

We consider an exogenons and deterministic investment cost, which depends on 

lifetime, T, and lead time, z/, among other factors. These two parameters, both considered 
deterministic, also affect directly the net present valne, ip, of a single investment. Moreover, 
we consider the disconnt rate, r, nsed by the hrm to also be exogenons and constant. In 
the backgronnd, we hx a complete probability space P), eqnipped with a hltration 

F = (-Fi)t>o- Let (Xt)t>o be an J^^-adapted ontpnt price process. An investment generates a 
random cash how-process of the form (/ (Xt))i>o, where / is a fnnction known to the hrm. 

The expected disconnted stream of fntnre cash hows, minns the initial investment cost, 
yields the net present valne 

pu-\-T 

ip{x]T,u) =—I{T,u) + / E {/} df, 0 < x < cxo, (1) 

J U 

where we have introdnced the conditional notation = {Xt|Xo = x}. The expression in 
(fTj) helps clarify the role of the lead time v as the time between the expenditnre I and start 
of operation and access to the cash how / (X^). 

The limitation on the amonnt of capacity that can be active at any given point in time, as 
implied by the assnmption of the role as a price taker discnssed above, natnrally introdnces 
the notion of a refraction time. That is, two consecntive investments have to be separated 
with at least a time T, the lifetime of the capital. With this being the only restriction on 
the investment strategy of the hrm, the valne v^^\x;T,u) of k consecntive investments can 
be formnlated as an optimal mnltiple stopping problem 


v^^\x-,T, u) 


snp E 




2=1 



0 < X < oo. 


( 2 ) 


The set is the set of all refracted stopping times r = (ti, T 2 , ■ ■ ■, Tk), i.e. Tj — Tj-i > T for 
i = 2,3 ..., k. Under this reqnirement, nothing prevents the hrm from having two consecntive 
investments operating with a seamless transition. This is becanse the decision to invest in 
fntnre capital can be made before the cnrrent investment expires. Provided that a least npper 
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bound in (Ej) actually exists, we can include in the stopping vector where Tp = Tp+i = • • • = 
Tk = oo, p < k. For such a stopping rule, can be interpreted as the value of the contract 
when not every exercise is being called. 

The tacit assumption that the supremum in ([2]) exists implies that T, z/) is in- 

tegrable and u)} = 0 for every choice of T and u. For brevity we will 

suppress the parameters T and u in '0, and /, unless such a dependence is specihcally 
required. Since the hrm is free to choose the timing of every investment, up to the condition 
of a refraction time T, the following proposition highlights the optionality embedded in the 
investment decision. A proof is provided in Appendix A 


Proposition 1. The value function k > 1, satisfies, 


v^^\x) = sup E 

tG5'= 


5^£!-’-"v+(Ay*) 


(3) 


2 = 1 


with = max{0, 

Proposition [1] reveals that it is never optimal to make the investment in the negative region 
of fiix). This is because the hrm has the hexibility to delay investment, and therefore, can 
always avoid value destruction. 

We now consider the formulation in ([T]) and (|2]) under the geometric Brownian motion 
(GBM) model. Specihcally, we model an output price process by the stochastic diherential 
equation (SDE): 

dXt = aXtdt + aXtdBt, Xo^(0,C)o), (4) 

where {Bt)t>o is a standard Brownian motion. The drift rate a and the variance parameter a 
are both assumed constant. We denote F = (Xt)t>o to be the hltration generated by B. 

From the hrm’s perspective, we assume that the investment cash how is the diherence 
between an uncertain output price W, modeled by (jl]), and a constant operational cost c. (In 
Section H] we show that such a formulation is equivalent to a cash how being the diherence 
of two GBMs, possibly correlated, through a change of measure.) Moreover, the hrm has the 
operational flexibility to temporarily suspend production if cost exceeds output price. Hence, 
the ehective investment cash how f{X^’^) at any future time t is given by 

= (5) 


can be 


where {X^ — cy = max{0,Xt — c}. As such, the expectation of f{Xy^) in 
seen as the price of a Europ ean call option on X with a strike price c and maturity t (see 
McDonald and Siegel f 1985h h namely. 


E I - c) + } = a;<F(d+(f))e“* - c<F(d_(f)), 
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where <h is the standard normal cumulative distribution function, and where 


d±{t) = 


X 


In ( — ) + a ± -a^ t 


/aVi. 


Applying this to ([1]), the reward function 'ijj{x) becomes 

ri'+T 


i>{x) = -I + 


(a:$(d+(f))e"* - c<h(d_(t))) e-^^dt . 


Under this setting, we recast the optimal multiple stopping problem as a sequence of 
optimal single stopping problems. More precisely, we express the value as 

v^^\x) = sup E • (6) 

The set S is the set of all F-stopping times and 

(x) =i:{x) + e-'^E = 0 . 

By construction, the admissible stopping times (ri,..., r^) are again refracted by at least T 
(years). By standard optimal stopping theory, the Snell envelope is con¬ 

structed to be the smallest supermartingale dominating the process 

ev ery k. For further details re garding this approach under a more general framework, we refer 
to ICarmona and Touzil (120081) . 


3. Analytical Results 

In this section we provide an analytical study for the optimal stopping problem in (|6]). Our 
main result is the characterization of the value function and optimal stopping rule, for every 
k > 1 opportunities to invest (see Theorem [1]). This leads us to develop an iterative approach 
to hnding the value functions n^^^(x) and the optimal exercise boundaries x^ (see Corollary 
[1]), which will be the foundation for our numerical algorithm to be presented in Section l4Tl 

In the context of real investments, let us discuss some conditions on the general reward 
function -0. First, we assume that 0 is continuous, increasing and sufficiently smooth. Ad¬ 
ditionally, we assume that there is a unique break-even point, xo > 0, where 'ip{xo) = 0, and 
0(a;) < 0 for x < xq and 0(a;) > 0 for x > xq. Finally, we note that 0(a;) is bounded by 
some affine function of x. Such a bound, together with an assumption that the discount rate 
r exceeds the drift rate a of the underlying process, ensures that perpetual waiting will lead 
to zero expected reward, i.e. linit^oo E {e“^*0(Xt)} = 0. Under these conditions, we will 
hrst analyze the real option problem with a single investment opportunity, which will be the 
building block for solving the optimal multiple stopping problem. 
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3.1. Optimal Single Stopping Problem 

We first consider the optimal single stopping problem 

= sup E . (7) 

TiG<S 


This problem is similar to the pricing of a perpetual American option with the payoff 'ip- 
Consider a candidate stopping time for problem ((Tj): 

Tx* = inf{f > 0 I > 0 }, 


with thr eshold xt > 0 . By the well-known Laplace transform of the hrst passage time of X, 
see e.g. fIShrevel . 12004 p.346), we have for x < 


where 





( 8 ) 



Note that the condition r > a implies that 7 > 1. If x > then = 0 and = ip{x). 

From (IH]) we see that a necessary condition for x\ to be an optimal stopping boundary is that 
x\ maximizes 'ijj{x)/x'^. The linear bound on 'ip{x) ensures the existence of such a maximum. 
The hrst order condition for a maximum a.t x = xl is given by 


d 'ip{x) xlpj^xl) — j'lpixl) 


(a:t)T'+i 


= 0 


Apj^xl) = 0 , 


(9) 


where we have introduced the operator notation A = (7 — x-^). The second order condition, 
sufficient to prove a maximum together with the hrst order condition above, can then be 
stated as 




X] 


dx'^ x^ 


d 

dx 


Ai/j{xl) 




< 0 


A 

dx 


Apj^xl) > 0 . 


Note that a maximum of 'ip{x)/x'^ at x\ is not a sufficient condition for an optimal stopping 
rule Txi for a general reward function ip. One would also have to prove that (Xt)) 

satishes the supermartingale property for this choice of The following lemma gives the 
sufficient conditions for Tx* to be an optimal stopping rule for problem (iTl). 

Lemma 1. Let ip : M’*' —)■ M 6e a reward function in the single stopping problem If xl is 
a global maximum for ip{x)/x'^ on M’*' and if 




(7 — l)ip'{x) — xip"{x) >0, X > xl, 


( 10 ) 
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then 


= ijj{x V x\) 

where v^^Hx) is continuous on M’*'. 


. I ^ 

lA ( - 
xl 


X e 


Remark 1. The first order condition A.'ip{x\) = 0 in ([9]), together with the condition that 
{d/dx)A.'ilj{x) > 0 , a: > bonnds the second derivative (convexity) of the reward fnnction 
'ip{x) for large x. Up to the condition of a maximnm of 'ijj{x)/x'^ at x = x^, the behavior of 
the fnnction '0(x) to the left of x = x^ is irrelevant. The conditions in Lemma [T] are therefore 
less restrictiv e than imposing that th e drift term of be monotone for all Xt, as 

presented in ( Dixit and Pind^^ . 1994 . pp. 128-130) as a part of the sufficient conditions for a 
connected stopping boundary at x = x^ for a perpetual American call on 'il^{x). 

Proof: From ([ 8 ]) we know that v{x), where 


n X = 


'ijj{xl) 


X < xl, 

X > X*, 


( 11 ) 


is a candidate for the solution. From the conditions on 'ip{x), it follows that, for x > x^, 


A'lp^x) = 

Tav,(x) = 

Looking at the drift ternj^] of e~''^v{Xt), 


7 ' 0 (x) — x'lp'lx) > 0 , 

(7 — l)' 0 '(x) — xip"{x) > 0 . 


( 12 ) 

(13) 


¥.{de-'^^v{Xt)] = e 


(xt)T' L 


+e 


—rt 


-r 07 -h -o'^7(7 - 1) 


1 


^{Xt<xi]dt 


-r^(Xt) + aXtij'iXt) + -a^X^nXt) 


^{Xt>xi}dt, (14) 


we observe that the hrst term vanishes identically, and the second term is non-positive fol¬ 
lowing from flT^ and flT^ . Hence, {e~'^^v{Xt))^^Q is a supermartingale, which implies that 

v{x) = E )} > E , t £ S. (15) 

The linear bound on ifj[x) implies that v{X^I^^) is integrable. Also, taking the limit 

t — 7 - cxD in (fT5|) and maximizing over r yields: 

v{x) > supE{e-""V’(^°’'')} • (16) 

T^S 


^For a general underlying process X, a localization procedure may be applied to eliminate the martingale 
term of de~^^v{Xt) upon expectation. 
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Conversely, choosing the specihc stopping time r = the process 

1 y \ 

is a martingale by construction and therefore 

v{x) = 

= E II < sup E } . 

J res 


t>0 


(17) 


The expressions in flT^ and ffT7|) together give the desired result, 

v{x) = v{x) = supE • n 

rG5 


3.2. Optimal Multiple Stopping Problem 

With Lemma [Hand Proposition [1] we can now state the main result. See Appendix B for 
a proof. 


Theorem 1. Let : M’*' —?■ M 6e a reward function with a break-even point xq. If A'ijj{x) is 
convex for X G (a:o,oo), with Kfj^x) increasing for large x, then, for every k>l, there exists 
an x\ > Xq such that 


v^’^\x) = V xl) 


where 


[ia(4 

)1 

, k>l, 

(18) 

L Wk 

/ J 



‘)(xp)}. 

(19) 


Moreover, the seguence {xDj^^^ is strictly decreasing, and is a strictly increasing 

seguence of continuous functions on M"*". Also, for any bounded subset D C M"*" there exists a 
constant K^, such that v^^\x) < K^, for x E D and k > 1. 


From a computational perspective it is convenient to introduce the auxiliary function 
through 

= AV>("na^)l{x>x*}. (20) 

With the conditions imposed on 'ijj{x) (bounded by a linear function, convexity of Aiplx) on 
(xo, oo) and A'0(x) being increasing for large x) one can, in a similar way as for xin the 
proof of Theorem m show that is an increasing sequence of continuous functions 

bounded on every bounded subset D C Given the function u^^\x) we can reconstruct 
the value function v^^\x) through 


v^’^\x) = x^ 


f 

V KP 


y ^ ^u^^\y)dy 


( 21 ) 


Note that since u^^i{x) = 0 for x < x^ there are no convergence issues in fl2T]) . As a conse¬ 
quence of Theorem [T] we have the following result. 
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(22) 


Corollary 1. The functions u^^\x), for k = 1,2,..., satisfy 

The boundary point x\ is the unique solution to 


i {Ai,{x) + e-m L„: > 0 . 


(23) 


Corollary [H together with equation fl2T]) . outlines the inductive algorithm used to hud the 
solution to ([6]). Details of the implementation are found in the next section. 

It now remains to consider the behavior of the solution to the optimal multiple stopping 
problem in the limit of inhnitely many exercise rights. That is, we will investigate the value 
function 




x) 


sup E 

TeS°° 




(24) 


Since (e where ax > ax — b > fj^x), is a supermartingale for all Xt and since the 

refracted stopping times (r„)„>i satisfy Tn > {n — 1)T, we have 




x^ 






} = 


ax 


n=l 


i=0 


1 — ■ 


(26) 


That is, v^°°\x) is bounded on every bounded subset of M'*'. Dehned through an integral in 
(1^ . this bound ensures continuity of v^°°\x) on every bounded subset of M’*'. Finally, we 
state the following convergence result. 

Proposition 2. The sequence v^^\x) converges uniformly to v^°^'^{x) on every bounded inter¬ 
val o/M+. 

Proof: Without loss of substantial generality, we assume that x < M. Let {Tf^)n>i be an 
optimal stopping rule for the value function v^°°\x) in fl2T|) . Employing a similar argument 
as in (j25j), with f{x) = ax bounding 'if{x), we obtain 




L n=l 

{ oo 

r 

< v^'‘\x)aM z 


—VT 

g '' 71,00 


^n=k+l 

^-{r-a)kT 


n=k-\-l 


1 — e“h-“)'r ■ 


(26) 


In the second step we have used the fact that is an admissible, but not necessarily 
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optimal, stopping rule for v^^\x). With x < M it follows from Theorem [Uthat is 

strictly increasing and bounded, and hence convergent on [0, M]. The uniform convergence 
of on X G [0, M] then follows from ([26]). □ 


4. Application to Infrastructure Investments 

In this section, we apply our analytical results in Section |3]to the infrastructure investment 
problem and discuss a numerical implementation. Moreover, we study the sensitivity of the 
result with respect to the key parameters of lifetime T and lead time as well as the process 
parameters. Lastly, the proposed framework is employed to compare investment scenarios 
with relatively short lifetimes and lead times to a scenario with both a long lifetime and lead 
time. Such a comparison will reveal a critical investment cost of the short-lived scenario below 
which it will be competitive with the long-lived counterpart. 

Recall from Section [2] the cash flow /(W), in which we have incorporated the flexibility 
to temporarily suspend production to avoid a negative cash flow, namely, 

f(X,) = (X, - c)+ . 


The reward function associated with this cash flow is 


'll}{x) = -I + 




(27) 


where 

/aVi. (28) 

Since /(W) is bounded by Xt, one realizes that there is a linear function bounding i/^(x). 
Investigating the derivatives of ^(x) yields 


d±(t) = 




^'(x) 



(x<I)(d+(t))e“‘ 


- c^(d_(t))) e-’'*dt 





(29) 


In the second step, we have identihed the derivative inside the integral as the Delta of the 
European call option. From fl271) - fl2^ we see that 'i/j{x) is continuous and increasing in x. 
Furthermore, since lima-^o'0(3^) = ~d < 0 there exists a unique break-even point Xq. Thus, 
'ijj{x) satishes all the conditions in the dehnition of a reward function presented in Section [3Tl 
By differentiation, we hnd that 

^[AV^(x)] = ^ ((y - 1)^^ + 


where (/>(x) is the density function of the normal distribution. It can be seen that A.'ip(x) is 
convex for all x > x', where x' = cexp [— (a + — |cr^) u\. In turn, if xq > x', which is 

ensured by imposing a lower bound on I, then XiIj(x) is convex on (xq, cxo), and consequently. 
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all the conditions for Theorem [T] are satished and the algorithm ontlined in Corollary [T] can 
be applied. 

If the operational flexibility to suspend operation is removed then the expected cash flow 
can be readily integrated and the reward function 'ipo{x) in this case would be 


ru+T 


=—I + / e E { — c) } df = a(i^, T)a; + T), 


where a{u,T) and b{u,T) are functions of the lead-time and lifetime only. That is, without 
the operational flexibility the reward function 'ipo{x) is clearly linear in x, rendering A'i/jo{x) 
convex on M"'", and thus Theorem 1 is also applicable to this case. 

4.I. Numerical Implementation 

The expressions for u^^\x) and v^^\x) in fl20|) and fl2T]) above form the basis of our nu¬ 
merical algorithm, whereby u^^\x), and x^ are computed iteratively for fc = 1, 2, 3,... 

The calculations below were carried out on a grid of 500 points regularly spaced between 0 
and Xraax, where the latter was determined by the process and model paramenters. The com¬ 
putation of the expectation in fl2^ is complicated by the fact that is not bounded 

on R"*". However, since our choice of reward function, '0(a:) in fl271) . is approximately affine for 
large x, so are A'lp^x) and Rather than truncating the distribution we can instead 

linearly extrapolate u^’^~^'>{x) outside the given grid. Since (x’Dk^i is decreasing, the necessary 
range of the grid depends mainly on x]". Specihcally, was chosen to be the upper bound 
on a two-sided 99.9% conhdence interval around x*. 


Xr, 


= exp 


In(xt) + [a--a^ 


T + 3.29(tVt^ 


This choice of Xmax, together with the number of gridpoints, was an acceptable compromise 
between having a large enough grid to ensure a linear behavior of without undue 

computational complexity. When comparing different scenarios, i.e. different values for T, u 
and I, the largest grid was used throughout. 

Using a trapezoidal method, the expectation in fl2^ was calculated as 


/ •^max / 

E x,T, a, a)dz + / {kz + m)g{z-,x,T,a,a)dz, (30) 


where kz + m is the linear extrapolation of for x > x^ax and where g{z]x,T,a,a) is 

the density function of the lognormal distribution. The upper limit x in fl30l) was chosen so 
that the support of the distribution contained a two-sided conhdence interval around x, for 
every x < Xmax- 

The boundary points xl were found using a simple bisection method on the convex function 
in (123|) . The calculations of u^^\x)^ and therefore also of x], and v^^\x)^ were terminated once 
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a tolerance e, defined by 


1)1 


e = 




(31) 


had reached £ < 10 “^. Such a tolerance yields a solution v^^\x) that is within 0 . 1 % of 
an accuracy that is likely good enough considering reasonable errors in the estimation of the 
underlying parameters. In the results below we denote by x*^ and the stopping 

boundary and the value function respectively at the termination of the algorithm according 
to the tolerance in fIM]) . 

Figured] demonstrates the convergence of the algorithm for the value function v^^\x) and 
the stopping boundary x\ respectively. Default parameter values used in the calculations 
below are given in Tabled! The choice of hxed cost parameters (/ and c), or their ratio J/c = 
10 , stems from observations in pertinent industries. Recent studies on energy investments and 
commodities report price volatilities in the neighborhood of 20% though they can occasionally 


experi ence short-term drastic increases (see lWestner and Madlenerl (120121 ) : lGeman and Ghana 


( 20091 11. Therefore, setting a = 20% corresponds to a relatively conservative choice as lower 
volatility favors the status quo. However, we will also examine a relatively high volatility 
(a = 40%) scenario in Figure[5l Additionally, the use of an effective discou nt rate [r — a] of 5% 
i s also in line with examples found in the literature in related helds, such as Fraver and Uludei^ 
( 200 il l. Overall, the choices of parameter values are, albeit reasonable, mainly for illustrative 


purposes. 

According to Theorem dl 


[v^^\x)) 


k>l 


is strictly increasing and the sequence of stopping 


boundaries is strictly decreasing. In Figure HI (left), we see that the value function 

increases monotonically with each iteration. After 50 iterations the tolerance e, dehned above, 
was less than 10“^. Moreover, we notice that the value function appears linear for large x. In 
Figured! (right), the stopping boundary x^. decreases rapidly from 0.85 to 0.44 and we mention 
that the break-even point was xq = 0.33. The convergence is clear even after 20 iterations. In 
particular, the stopping boundary x* = 0.85 (see (|9|)) from the hrst iteration helps us dehne 
the upper bound x^ax of the grid. 

The stopping boundary x\ is the price level at or above which the hrst investment should 
be made, given the option to make k — 1 more investments later, each separated in time by at 
least the lifetime T. Since hrms rarely have an imposed number of investment renewals, the 
boundary x^ is of primary interest. Contrasting the value of x^ to x^ reveals the impact of 
including future investment options on current investment decision. In Figured! the number 
of exercises (iterations) is A; = 50, which together with the lifetime of T = 5 years, gives a time 
horizon of more than 250 years. Such a time horizon is practically inhnite under reasonable 
circumstances. On the other hand, with a large /c, the number of remaining investment 
opportunities should have a smaller impact on the hrst investment timing. This is evidenced 
by the convergence of x\ to a constant level as k increases. 

From our numerical tests, we hnd that the refraction time (lifetime) T strongly inhuences 
the speed of convergence. This is intuitive due to the discount factor in the dehnition 
of u^^\x) in (l2^ . reducing the differences over iterations. On the other hand, the lead time 
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Description 

Parameter 

Value 

Lifetime 

T 

5 

Lead time 

V 

1 

Investment cost 

I 

1 

Operational cost 

c 

0.1 

Discount rate 

r 

10% 

Drift rate 

a 

5% 

Volatility 

a 

20% 


Table 1: Default parameter values used in the calculations. 


V has much less bearing on the rate of convergence since it only affects the hrst investment 
timing. We will further examine the impacts of lifetime and lead time in the next subsection. 




Figure 1: (Left) Convergence of the value function, for A: = 1,..., 50. (Right) The stopping 

boundary x^. decreases rapidly over iterations. 


4-2. Sensitivity Analysis 

We will hrst investigate the sensitivity of with respect to the process parameters a 
and a. Later we investigate the sensitivity of both and the value function with 

respect to the main model parameters of lifetime T and lead time u. 

Low values of the drift rate a corresponds to a higher ‘effective’ discount rate r — a, some¬ 
times called the convenience yield. This decreases the present value of future investment, 
which explains the minor difference between and for small a, see Figure [2] (left). Con¬ 
versely, a higher drift rate of the underlying enhances the value of future investments. This 
drives a greater wedge between the stopping boundary xl of a single investment compared 
to the stopping boundary of multiple consecutive investments for large a. Consequently, 
this emphasizes the importance of including future investments in current decisions in envi¬ 
ronments with a high drift rate. Note that the same result is to be expected if decreasing the 
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discount rate r, still with the condition that r — a > 0. In the figure, we observe that as a 
approaches r (10%) the exercise boundaries /c = 1, 2, 3, increases rapidly. This is intuitive 
because theoretically the optimal exercise boundary would be infinite for a > r. The same 
phenomenon also occurs for when a is very close to r. 




Figure 2: Sensitivity of the optimal stopping boundaries w.r.t. the drift rate a (Left), and w.r.t. 
the volatility a (Right). 

In Figure [2] (right), we see that the optimal exercise boundary increases with volatility 
cr. This suggests that in a more volatile environment the firm will demand a higher output 
price level in order to enter the market. We observe that the increasing pattern holds for 
finite /c, k = 1,2,3, as well as for the infinite case. 

Now, we turn to examine the impact of lifetime T, the parameter of main interest, on the 
stopping boundaries xl and value v^^^(x), with the same fixed cost 1 = 1. First, in Figure |3l 
we observe a minor difference between x\ and x*^ for long lifetimes (for T > 20). Intuitively, 
the incremental value of an additional investment 20 years or more from now is minimal due to 
the discount factor in fl2^ and (|23|) . Therefore, for very long lifetimes, future investment 
decisions will not significantly influence the current decision to invest. However, for very short 
lifetimes, T < 2, there is a substantial difference between the optimal exercise levels with one 
and infinite investment opportunities (i.e. x\ vs. x*^ for small T). In an intermediate regime, 
with lifetimes of 5 to 15 years, we observe that including only one or two future investment 
decisions will significantly affect the decision regarding the first investment. In both finite 
and infinite cases, the optimal exercise boundary decays rapidly with respect to lifetime. For 
instance, x*^ is almost fiat for lifetime of 10 years or longer. 

Considering the value v^°°\x) of the multiple investment scenario for different lifetimes 
we observe a trend similar to the one regarding stopping boundary. The marginal impact of 
adding one year of life to short-lived capital is significant. However, this impact drastically 
decreases for longer-lived capital. For instance, as seen in Figure |3l the value of multiple 
investments in capital with a 25 year lifespan is virtually the same as capital with a life of 5 
years, at the same investment cost I and the same lead time u. Even though the investment 
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Figure 3: (Left) Optimal stopping boundaries decrease as lifetime T increases. (Right) The value 
v^°°\x,T) is increasing with respect to the underlying price levels x and lifetimes T. 


cost is the same, the increased flexibility in timing future investments in the scenario with the 
shorter lifetime makes them almost equally attractive. 


1.5 



m 


Ql-.-.-.-^-1 

0 1 2 3 4 5 

Lead-time, z/ 



Figure 4: (Left) Optimal stopping boundaries for different lead times v with the same fixed cost 
1 = 1. (Right) The value v^°°\xo,Ti,u) evaluated for different lead times u at xq = 0.5, and for 
Ti = 5,10,15,20. 

As seen in Figure H] (left), the lead time does not substantially influence nor the 
difference between x\ and x*^ at fixed v. In other words, varying lead times (within the 
given range) will barely affect the firm’s decision to invest. This is expected since seamless 
consecutive investments are possible regardless of the length of lead time. On the other hand, 
a longer lead time delays revenue generation relative to the outlay of the investment cost 
I. Increasing lead time therefore decreases the net present value of every future investment. 
This explains the decreasing trend of the value function with respect to lead time in 
Figure H] (right). In summary, the firm will realize a lower value with longer lead times of 
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the investments, even thongh the corresponding exercise bonndary changes only marginally. 
Displaying the value at different lead times and at a hxed price level but for several 
values of the lifetime T again shows the diminishing returns of adding lifetime to capital. 

From Section [2l we know that the reward function -0 can be interpreted as the sum (inte¬ 
gral) of European call options on the uncertain output price X, with strike c and maturities 
ranging over the lifetime. Since increasing the strike price of a European call decreases its 
value, the cost parameter c has the same effects on the value which in turn raises the 
stopping boundary x^. By similar reasoning, the same holds for I. 

Lastly, we examine the value of the option to temporarily suspend production for invest¬ 
ments with short and long lifetimes. Two scenarios are compared in Figure [5l short-lived 
investments (T = 5, u = 0.5) and long-lived investments (T = 25, = 5) at two levels of 

volatility a = 20% and 40% respectively. Note that in both cases an investment cost of / = 1 
was used, which explains a relatively higher stopping boundary for the short-lived capital. In 
this example, the operational flexibility appears to barely affect the investment decision in the 
short-lived scenario regardless of volatility. For low volatility the same indifference appears 
in the long-lived scenario as well. Figure |5] (right). However, in contrast to the short-lived 
scenario, the operational option leads to earlier investment timing under high volatility in the 
long-lived scenario. Comparing the left and right panels of Figure [H we see that a shorter 
lifespan yields higher thresholds for future investments. 




Iteration, k 


Figure 5: (Left) Stopping boundary for a short-lived scenario where capital has a lifetime T = 5 years 
and a lead-time n = 0.5 years. (Right) Stopping boundary for a long-lived scenario where capital 
has a lifetime T = 25 years and lead-time n = 5 years. 

This operational flexibility feature is typically associated with competitive and unregulated 
markets, whereas many of the relevant industries for this paper traditionally have operated 
in regulated environments. Moreover, it could be argued that the conventional ‘bigger-is- 
better’ paradigm has historically invited regulations on the market since monopolies are often 
observed in industries with a large production scale. However, it is conceivable that these 
conditions may change as new technological advances allows for more competitive distributed 
small-scale production models. 
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As an example, these changes are occurring and subsequently making operational flexibility 
available in the electric power generation. In the United States, more than half of the states 
have moved towards deregulated power markets. While power producers face many different 
costs the dispatch order of the various available generators in a given market is to zeroth 
order a function of operational cost. That is, it is up to the individual producers to decide a 
price point where operation commences. It should be noted that the interpretation of a unit 
not operating is dependent on the time scales involved. For instance, hour-to-hour decisions 
typically implies putting generators into an synchronized reserve capacity called spinning 
reserves, whereas seasonal interruptions correspond to complete cold stop. Some technologies, 
like nuclear power, currently have a relatively low operational cost and are thus ahead in the 
dispatch order to meet demand (also called base load generation). In general, the value of 
the operational flexibility option is expected to be higher if the production involves lowers 
costs of suspension and resumption. Lastly, we remark that our standing assumption that the 
Arm acts as a price taker would have to be revisited if the resulting change in supply should 
significantly impact the market price. 

4-3. Critical Investment Cost 

Many, if not all, of the fundamental process industries, such as e.g. energy, water and 
petrochemicals, have followed the trend of ‘bigger-is-better’ over the past century. Experience 
from a given industry gives us the investment cost (per unit of capacity) hargei as well as the 
lifetime Tiarge and lead time viarge of large-scale capital in the current paradigm. Transitioning 
to a paradigm of small-scale, mass-produced, and modular equipment will give rise to a new 
parameter ensemble, {Ismaii,Tsmaii, ^smaii}, where we a priori only infer shorter lifetimes and 
lead times. By considering an infinite time horizon, we can use the framework in this paper 
to And the critical investment cost Icrit that would render a small-scale approach competitive. 
That is, for a given reference ensemble {harge, Tiarge, Carge}, we And for every choice of Tsmaii 

and Usmaii an Cru such that vl'^hi > provided that Cmaii < Icrit- 

The observed historical trend of increasing unit sizes seemingly suggests that total cost 
decrease with size. Nevertheless, we will here let the operational cost c be independent of the 
parameters I, T and u. To motivate this setting for our examples below, we consider the total 
capacity to be installed at a single location, so ancillary costs, such as transportation costs 
of inputs and outputs to and from the plant, administrative costs, cost of security etc., ought 
not depend on the granularity of the hardware inside this black box. This leaves primarily 
two potentially size-dependent factors that can affect operational costs; labor and efficiency. 
Historically, the necessary amount of operational labor increases with the number of individual 
units employed. Operating fewer and larger units therefore increases labor productivity and 
lowers labor cost per unit output. However, advances in automation make it possible today to 
decouple the amount of required labor from the number of individual units at reasonable cost. 
As to efficiency, from a physical perspective one can in some cases claim that larger units have 
lower dissipative losses (e.g. through friction and unwanted heat transfer) and hence are more 
efficient and therefore less costly to operate. On the other hand, it is not always clear how 
important these effects are. For instance, a statistical analysis on the size-dependency of cost 
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in four different electricity generating technologies showed that o nce labor cost is remo ved. 


unit size is not a signihcant factor affecting total operational cost fjPahlgren et al.l. 120131) . 


As an illustration, we can compare a single-cycle thermal power plant to an internal 
combustion engine, both performing fundamentally the same task of converting chemical 
energy into mechanical work. Under reasonable circumstances, they can do so at comparable 
efficiencies. The car engine is mass-produced on the order of days. The power plant, on 
the other hand, is typically not ready for operation until several years have passed since the 
decision to invest was made. Moreover, the power plant is designed to last for decades while 
the car engine presumably will have a lifetime on the order of years under constant operation. 
How much would one be willing to pay for an engine that is fully automated, retrohtted to 
run on the same fuel and equipped with a generator to produce electricity? Such a mini-power 
plant can reasonably be assumed to incur similar levels of operational costs per kWh produced 
as its large-scale counterpart. This leaves investment cost, lifetime and lead time as the main 
distinguishing features from the large-scale power plant. 

We start by analyzing the reward function in fl27|) for large x. It can be seen that $((i±) ~ 1 
for X ^ c, and therefore, is asymptotically affine with 


»i/+r 


'ip{x) —1+ — c) e '^^dt 

J U 

^-(r-a)u 


a 


-(1 
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)c ) = ax — b. 


(32) 


For large enough values of x the function is affine as well and 

k 

- l)aa: - , 


u'- '[X] = 


i=0 


where E = Using (l3^ to substitute for a, we obtain for large x\ 


u^°^\x) 


3 —(r— 


-X 


■yh 


r — a 


1 - ’ 


where the slope is independent of the lifetime T and decreasing in the lead time u. From the 
dehnition of in (12(11) it follows that the same properties can be ascribed to the value function 
yd^') (^x; T, u) for large a: S> c. Consequently, comparing different scenarios, i.e. different T, u 
and /, the scenario with the shortest lead time will always have the higher value for large 
enough values of x. 

Now we compare the value of small-scale capital, here characterized by Tsmaii < 5 years, 
k'smaii < 3 years, against a benchmark of Tiarge = 25 years, Piarge = 5 years, representing 
traditional large-scale capital. Furthermore, as a point of reference we set harge = 1, which 
together with the constant operational cost c = 0.1 (assumed the same for every choice of /, 
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T and u) provides a relative monetary scale. 




Figure 6: (Left) The value of a large-scale investment scenario {Tiarge = 25 and Vsmaii = 5) together 
with small-scale analogues {Tsmaii = 3 and Ugmaii = 0.25) at different investment costs. The figure 
verifies the result that the scenario with the shorter lead time has the higher value for large prices, 
X. (Right) Displaying the ratio of the previous value functions clearly 

reveals that ^^(x) > ^very x. This suggests the existence of a critical investment 

cost Icrit such that the value of a small-scale investment scenario exceeds the traditional large-scale 
counterpart for any price x, as long as Ismail < Icrit- 

In Figure [6l the value of the large-scale parameter ensemble is displayed along¬ 
side with Tsmaii = 3, lygmaii = 0.25 for three different investment costs, Ismaiij = 

0.5,1,1.5. With their shorter lead time, the values j(^) exceed v\^g^{x) for 

large values of x, verifying the analysis above. 

From Theorem [T] we know that v^'^\x) is proportional to x'^ for small enough values of x. 
This explains the constant appearance of the ratio Vg^lii j{x)/vl^g^{x) in Figure [6]for small x. 

The value of the scenario with the lowest investment cost clearly exceeds the value of 

the ‘large’ scenario for all prices x. Indeed, given Tiarge J^iarge and Iiarge we can find a critical 
value IcritiTsmaii, i^smaii), such that the value of the ‘small’ scenario is greater at any price 
level, as long as Ismaii < Icrit- This leads to the contour plot in Figure [71 which displays the 
critical values for different lifetimes, Tsmaii and lead times, vsmall- In Section l4)2l it has been 
observed that the value increases with lifetime (see Figure |3]) and decreases with lead 
time (see Figured]). This helps explain the trend of Icrit with respect to T and u in Figure 
[71 Precisely, in this domain of short lifetimes, a reduction in lead time yields a higher critical 
investment cost, which in turn increases the competitiveness of the small-scale approach. 

As an example, with Tsmaii = 2.5 years and Vsmaii = 0.3 we can infer from Figure [7] that 
Icrit = 0.5. That is, despite a difference of a factor 10 in lifetime {Tiarge = 25 years), the 
investment cost is required to differ only by a factor of 2 between the short and long lifetime 
scenarios in order for the short-lived one to be preferable. 
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Figure 7: Critical investment cost Icrit{T, u) of a single small-scale investment, with the given T and 
u, in order for the value of multiple consecutive such investments to exceed the value of multiple 
consecutive large-scale investments at any price level. Since the value function, for any choice of 
T and is decreasing in the investment cost, we have > v^°°\liarge,Tiarge,i^iarge)i 

provided that I < Icrit{T, v). The large-scale investment was characterized by the parameter values: 
harge = 1 , Tiarge = 25 years and i^iarge = 5 years. 


The closest resemblance of an estimation of Icru using a simplistic NPV argument without 
any optionality would be to compare the discounted investment costs over an inhnite horizon. 
That is, we can assume that capital is replaced every Tiarge (resp. Tsmaii) years under the long 
(resp. short) lifetime scenarios. For simplicity, disregarding features like lead times, opera¬ 
tional flexibility, stochastic prices, and optimal multiple stopping, we equate the discounted 
costs over an inhnite horizon to get 

OO OO jf 

E ThTio^j^g^ j — ^ TkTg^Q^ii jf crit 

^ J-large — / ^ ^ J-crit ^ j 

h=n h=n 


_ g tTsjykiII 


where denotes the critical investment cost of the short-lived scenario under the NPV 
argument. Using the same example as above, with Tsmaii = 2.5 years, Tiarge = 25 years, and 
harge = 1, we hud that = 0.24. Comparing to the critical investment cost Icrit = 0.50 
above, the value of the optionality in this particular example permits twice the investment 
cost that would be suggested by standard net present cost arguments. 

Lastly, we note that, in the example with the car engine and the power plant, the cost 
per kW of capaci ty of the car engine is a lmost two orders of magnitude less than that of 
the power plant f Larminie and Did3. 20031) . Clearly, the critical costs suggested in Figure 
[3 are not nearly as dramatic. This suggests great potential in abandoning the customized. 
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large-scale investments in favor of mass-prodnced and modular capital. 


5. Conclusions 


We have developed a framework for valuing repeated infrastructure investments and com¬ 
paring traditional capital with long lifetimes and lead times to a mass-produced version with 
a much shorter lifespan and lead times. We found that including future investments in present 
valuations signihcantly affects the exercise boundary, especially when the individual invest¬ 
ment is short-lived. Also, the marginal beneht of increasing lifetime rapidly diminishes when 
considering multiple investments. Furthermore, this framework allows for a determination of 
a critical investment cost of capital with lifetimes and lead times that deviate from indus¬ 
try standards. One of examples (see Figure!^ reveals that reducing lifetime of capital from 
a typical 25 years to 2.5 years need only be accompanied by a decrease of a factor of 2 in 
investment cost in order to be superior, in overall value terms. 

For our analysis and tractability, we have worked with a lognormal stochastic cash flow. 
Working with the lognormal process, we obtain analytical value function for the single stopping 
problem and analytical results like Theorem 1. Alternatively, for more realistic representation 
of commodity price dynamics, other processes, e.g. with mean reversion and jumps, can 
be used, though these settings may not be as amenable to mathematical analysis and the 
numerical implementation can be much more challenging and cumbersome. For instance, 
under the exponential Ornstein-Uhlenbeck price dynamics, the single stopping problem does 
not admit a closed-form solution, and t he op timal exercise threshold can only be found from 


an implicit equation fsee iLeung et al.l (120141) 1. However, under mean reversion, the value of 


an investment with considerable lead time and lifetime is expected to generate a cash flow 
close to the hnite long-run mean. Shorter lived, and more quickly deployed capital would, on 
the other hand, be more suited to both exploit positive deviations from the mean, and also 
avoid periods of low prices. In the mean reversion framework, we would therefore anticipate 
that the critical investment cost, IcHti of the small-scale investment to achieve parity between 
the small-scale and large-scale scenarios to be higher than those found in Figure [71 This 
suggests that a small-scale and modular approach may be attractive to large infrastructure 
investments. 

There are a number of directions for future research. First, one can incorporate the hrm’s 
av ersion to risk and ambi guity associated w ith the stochastic investment returns, as discussed 


m [Henderson] (120071) and I.Taimnngall (120111) in the case of single investment. Moreover, addi¬ 
tional information on the technological change, industry outlook, and future investment cost 
will enhance the decision analysis. For instance, intuition suggests that frequent technological 
advancement is more easily harnessed with a more rapid turnover since outdated technology 
can be abandoned without sacrihcing investments with a long remaining horizon. Similarly, 
regulatory risk is less of an issue with shorter investment windows. This issue is very much 
in focus in the energy sector today with the looming threat of climate change and related 
concerns about carbon emissions. For instance, burning cleaner than coal and oil, natural 
gas is touted as a bridge fuel to cleaner technologies in the future. However, in the current 
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paradigm of large-scale units with long lifetimes, the dramatic increase in natural gas-hred 
capacity in recent years means a commitment to a fossil fuel for at least another quarter of a 
century. 
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Appendix A. Proof of Proposition [T] 

For every hxed stopping rule r = G we dehne a subsequence s < A;, of 

by recording only those stopping times at which the reward is non-negative: 

{fj} = {n I > 0 ; f = 1,..., A:} . 

In the case of hnite k we can append the subsequence with k — s inhnite stopping 

times, fs+i = ■■■ = fk = oo to create the full sequence f = Since the stopping times 

within r are refracted with at least the constant time T, the same is true for those in f by 
construction, and therefore f G iS^. 

By avoiding those stopping times with a negative reward, the total discounted reward: 

k 

f e 5‘, 

is dominated by gk{x;f,'ip) in expectation, namely, 

E{gk{x;f,^p)} <E^^gk{x;T,'ijj)Y (A.l) 

In addition, when r is taken to be f, we have the equality gk{x;f,Y) = gk{.x;f,'il!~^), almost 
surely. That means that maximizing over the stopping rules r for the original problem v^^\x) 
will achieve the upper bound (RHS of ([3])) with the non-negative reward . In summary, 
we obtain 

v^^\x) = sup E{gk{x]f,ip)] = sup E igk{x;f,'ilj^)\ = sup E {gk{x;f,Y^)} . □ 

Tg.s'' fes* ^ ^ t£S'‘ 

Appendix B. Proof of Theorem [1] 

Proof: Be the dehnition of the reward function with a break-even point xq, we know that 
Yix) <0 and > 0 for x < Xq, and > 0. This implies that 

A'ljjix) = 7 ' 0 ( a ;) — x^'^x) < 0 , x < Xq. 
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Particularly, since h.'ip^x) is convex on [xq, oo) and increasing for large x, there is exactly one 
solution, xl, to A'ip{xl) = 0, and furthermore, {d/dx)K%lj[x\) > 0 for x > xi. Recall from 
Lemma [T] that 

v^^\x) =y xl) , (B.l) 

where In addition, 

^/;(2)(a;) = ^(1)(a:) + } . (B.2) 

Since AE{ 5 f(Xi)} = E{A 5 f(Xt)} for any integrable function g, we apply (IB.2P to get 




(B,3) 


where in the second step we have used the fact that An’^^)(x) vanishes in the continuation 
region of i.e. for x < x*. Since Aip^^^x) is assumed convex on (xq, cxo), the expectation 

E |A'0(^)(X°’'^)l{Xt>a;*}} is also convex on (xo,cxo). Being the sum of two convex functions, 
A'^(^i(x) is also convex on (xo,C)o). Moreover, since A'^(^i(x)l{ 3 ;> 3 ,*} is an increasing function 
(and strictly positive for x > x^), flB.3D implies that Aip^‘^\x) is increasing for large enough x. 

We observe from fIB.ip and flB.2D that il)^‘^\x) is a continuously differentiable increasing 
function with lim 3 ;^o = liiHx-i^o < 0. Furthermore, if 'iIj^^\x) is bounded by 

/(x) = ax, for some a > 0, then one can show that n*'^^(x) < ax, and therefore 

< ax + e“’’^E [aX^’^'^ < a (l + x. (B.4) 


This ensures the existence of a maximum at x = X 2 to the function '0^^^(x)/x'^, and conse¬ 
quently also the existence of a solution to Aip^‘^\x) = 0. Also, Proposition [1] implies that 
X 2 > Xo as it is never optimal to exercise with negative payoff. The convexity of Aij:^‘^\x) on 
(xo, 00 ) ensures that there is exactly one such maximum, uniquely dehned by 


A'iIj^‘^\x* 2 ) = 0, and -^A%jj^‘^\x) >0, x > x^. 
Hence, applying Lemma [U again we obtain 


n*'^^(x) = ijj^‘^\x V X 2 ) 



Also, from the bound on ip^‘^\x) in (IB.4P we infer that v^‘^\x) < a (l -1- e x. 

Repeating the argument above, we can similarly show the existence and uniqueness of a 
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stopping boundary xl, for every k > 1, and also derive an upper bound for each 

'k-l 


v^^\x) < ax I ^ ( 


, —(r— q) 2 T 


(B,6) 


, i=0 


Next, we proceed to prove by induction that which together with the pre¬ 

vious part of the proof would imply that x\j^-^ < x\. As previously remarked, A'0(^)(x)l{a,>3;*} 
is positive for x > x\. From fIB.Sp it follows that 

A'ip^‘^\x) - A-il:^^\x) = e“^^E > 0, 

Now, assume that AiIj'^^\x) > Aip^^~^\x). From the dehnition of ip^^\x) in flT^ we have 

(x) - AV^(^) {x) = e-'^E { Au(^) } 

aO,x\ -■ 


= e-'’^E 

> 0 . 

As for the monotonicity of the sequence we first show that Since 

> 0 from Lemma [H it follows from (IT^ that 


'ijj^'^\x) — 'ijj^^\x) = e ^^E > 0. 

By the fact that X 2 < x\ and il)^‘^\x) > the inequality 

■y(2) _ > Tjj^^\x) = 

holds for X > x\. Moreover, since ij:^‘^\x)/x'^ is maximized at x = X 2 , we see that 

(x"i) (Aj ie(i2,i;). 


(B,6) 


Similarly, for x < X 2 , 


X 


xl 


> i’‘-^\xV. 


X 


Xl 


> ^p^^\xl) 


X 


Xl 


= n^^^(x), X < X2 . (B.7) 


Hence, we have shown that gy induction, we obtain 

= e-"^E > 0. 

With this inequality we can follow the steps from (IB.61) to (IB.7P to arrive at (x) > (x). 


25 









Finally, the inequality in 
for any k: 


flB.5l) implies that the value function admits the following bound 
aM 


v^^\x) < 


1 — ’ 


X < M. 


□ 
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